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Abstract. Phase transitions in non-equilibrium steady states of 0(n)-symmetric models with re- 
versible mode couplings are studied using dynamic field theory and the renormalization group. 
The systems are driven out of equilibrium by dynamical anisotropy in the noise for the conserved 
quantities, i.e., by constraining their diffusive dynamics to be at different temperatures T" and T 1 - 
in d\\- and eij_-dimensional subspaces, respectively. In the case of the Sasvari-Schwabl-Szepfalusy 
(SSS) model for planar ferro- and isotropic antiferromagnets, we assume a dynamical anisotropy 
in the noise for the non-critical conserved quantities that are dynamically coupled to the non- 
conserved order parameter. We find the equilibrium fixed point (with isotropic noise) to be stable 
with respect to these non-equilibrium perturbations, and the familiar equilibrium exponents there- 
fore describe the asymptotic static and dynamic critical behavior. Novel critical features are only 
found in extreme limits, where the ratio of the effective noise temperatures T" /T 1 - is either zero 
or infinite. On the other hand, for model J for isotropic ferromagnets with a conserved order pa- 
rameter, the dynamical noise anisotropy induces effective long-range elastic forces, which lead to 
a softening only of the d± -dimensional sector in wavevector space with lower noise temperature 
T x < Tl The ensuing static and dynamic critical behavior is described by power laws of a hith- 
erto unidentified universality class, which, however, is not accessible by perturbational means for 
d\\ > 1. We obtain formal expressions for the novel critical exponents in a double expansion about 
the static and dynamic upper critical dimensions and du, i.e., about the equilibrium theory. 

PACS. 05.70.Ln Non-equilibrium thermodynamics, irreversible processes 64.60.Ak 
Renormalization-group, fractal, and percolation studies of phase transitions - 64.60.Ht Dy- 
namic critical phenomena 



1 Introduction be developed a bit further by noting that magnetostatics 

can be described in terms of interactions between currents 

The equilibrium properties of a generic system in contact and thus asking if NESS could also be described in terms 

with a heat bath are entirely determined by interactions of some effective interactions, characterizing the station- 

(the dynamics plays a role only in providing the necessary ary probability distribution. 

mixing in phase space). In contrast, dynamics is important The search for such effective interactions has been go- 
in non-equilibrium steady states (NESS) where competing ing on for some time. The simplest version of this search 
dynamics (i.e. contacts with more than one heat baths or consists of taking the logarithm of the steady-state dis- 
driving fields) generate fluxes of energy, mass, etc., or, tribution function for small-size systems and looking for 
equivalently, yield non-zero steady-state probability cur- the dominant interaction in the emerging effective hamil- 
rents in phase space. This difference between equilibrium tonian. Unfortunately, this approach did not turn out to 
and NESS has often been illustrated by using an electro- be very useful since, generically, all the interactions, of 
dynamical analogy jlj: equilibrium = electrostatics while any range, which are consistent with the symmetries of 
NESS = magnetostatics. The above analogy can actually the system are generated — and the lack of significant 
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differences in the magnitude of the couplings renders the 
identification of the dominant interactions impossible. 

A somewhat more sophisticated approach is based on 
the extention of universality concepts to phase transitions 
in NESS. The idea here is that the universality class of 
a phase transition in NESS provides information about 
the range of interactions generated by the competing dy- 
namics and, furthermore, also gives the exponent of the 
long-distance, power-law decay in case of long-range inter- 
actions. This approach has yielded some very interesting 
results in connection with non-equilibrium perturbations 
imposed on the relaxational models A (non-conserved or- 
der parameter dynamics) and B (conserved order param- 
eter dynamics) , using the terminology of Ref. Q . 

In particular, phase transitions in systems with model- 
A type of dynamics have been shown to be remarkably 
robust against the introduction of various competing dy- 
namics which are local and do not conserve the order pa- 
rameter ||, including competing dynamics which breaks 
the discrete symmetry of the system [Q . This robustness 
of the critical behavior also persists when the compet- 
ing dynamics comes from a reversible mode coupling to 
a non-critical conserved field [^) . Thus there is large class 
of non-equilibrium steady states where the character of 
interactions is not modified by the presence of thermody- 
namic fluxes or, in other words, the probability currents 
in phase space become irrelevant for the large-scale be- 
havior. Thinking in terms of the electrodynamic analogy, 
one may say that this corresponds to the magnetostatic 
problem of a steady electric current along a straight line, 
which can be reduced to an electrostatic problem after an 
appropriate coordinate transformation. 

Another class of competing dynamics is obtained when 
model-B type dynamics competes with external drive jl| 
or with local, anisotropic, order-parameter conserving pro- 
cesses H]t]J§. In these cases, long-range interactions do get 
generated in the NESS and, furthermore, a common fea- 
ture of these interactions is that their angular dependence 
resembles the form of elastic or uniaxial dipolar forces 
(and for some systems even the power-law decay with dis- 
tance is that occurring for dipolar interactions @,@,||). In 
terms of the electrodynamics analogy, one might say that 
the fluxes in the NESS of these systems are equivalent to 
loops of electric currents which interact, in the first-order 
multipole expansion, via (pseudo-) dipolar forces. (How- 
ever, when drawing this analogy, one should be aware that 
while successive terms in the electrodynamics multipole 
expansion become more and more suppressed, the effec- 
tive long-range forces appearing in the NESS of models 
with conserved order parameter are very dominant.) 

There are, of course, examples which do not fit eas- 
ily into this heuristic straight-line current and loop cur- 
rent classification, which therefore should not be taken 
too far; e.g., competition of the usual model A with non- 
local dynamics |9| or with linear coupling to a conserved 
field |1(J generates isotropic long-range interactions. Non- 
thermally driven steady states occuring, for example, in 
the prominent non-equilibrium universality classes of di- 
rected percolation flllfl and Kardar-Parisi-Zhang surface 



growth [ p2[ probably cannot be described in terms of such 
effective interactions at all. Furthermore, there are driven 
systems where the emerging effective interactions are ex- 
tremely long-ranged in the sense that the potential is non- 
integrable and, consequently, the system does not display 
thermodynamical behavior in the usual sense Jl3|,[l4| . 

Nevertheless, since a large variety of competing dy- 
namical processes have yielded a surprisingly small num- 
ber of effective universality classes, we feel that it is worth- 
while to continue the exploration of NESS through study- 
ing non-equilibrium phase transitions and thus deducing 
effective interactions. The investigation of the stability of 
the known equilibrium dynamic universality classes is im- 
portant also from another viewpoint. Namely, in many 
experiments probing critical dynamics, it is by no means 
a trivial issue to maintain thermodynamic equilibrium, 
as relaxation times become very long close to a second- 
order phase transition as a consequence of critical slow- 
ing down. For the interpretation of experimental data, it 
might therefore be important to know if the dynamical 
system is driven to a NESS that is characterized by scal- 
ing behavior which is distinct from its equilibrium coun- 
terpart. 

In the present work, we follow up a previous inves- 
tigation H of NESS generated in 0(n)-symmetric sys- 
tems subject to non-equilibrium perturbations. Here, ei- 
ther a non-conserved order parameter is dynamically cou- 
pled to non-critical conserved quantities — this defines 
the n-component Sasvari-Schwabl-Szepfalusy (SSS) model 
[ p^Ul^ |, which incorporates model E for planar ferromag- 
ncts (with n = 2) []l7| and model G for isotropic anti- 
ferromagnets (n — 3) |Q as special cases; or the or- 
der parameter and the generators of 0(n) are identical, 
which is realized for n = 3 in model J for (idealized) 
isotropic ferromagnets ]l9] . (For a review on more realis- 
tic dynamics for ferromagnetic systems that includes the 
effects of dipolar interactions, see Ref. [§0|.) This field- 
theoretic renormalization-group (RG) study found that 
spatially isotropic violation of detailed balance generically 
leaves the equilibrium fixed point stable, thus indicating 
that the steady-state fluxes involved here do not generate 
any long-range interactions. Since dynamical anisotropy 
appears to be an essential ingredient for the generation of 
pseudo-dipolar effective interactions in NESS, it is natural 
to ask if this was the case for the above 0(n)-symmetric 
systems as well. More precisely, we allow for dynamical 
spatially anisotropic noise for the conserved quantities. 
For the anisotropic non-equilibrium model J, this means 
that we allow for different effective temperatures T" > 
T , respectively, governing the order parameter noise in 

and d± -dimensional subspaces (dy + d± = d). Essen- 
tially, this means that there is a non-zero heat current 
flowing from the "hotter subspace" (d|| dimensions) into 
the <ij_-dimensional "cooler subspace" . For the anisotropic 
non-equilibrium SSS model, a similar distinction applies 
for the noise in the Langevin equation that describes the 
dynamics of the purely dynamically coupled conserved 
field, and we may then explore its influence, and that of 
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the ensuing effective heat current, on the non-conserved 
order parameter dynamics. 

In Sec. ^, we shall discuss our basic model equations for 
the spatially anisotropic non-equilibrium generalizations 
of Langevin dynamics appropriate for second-order phase 
transitions in 0(n)-symmetric systems. We start our in- 
vestigations with the anisotropic non-equilibrium version 
of the SSS model in Sec. ||, and shall find that for this sys- 
tem of coupled Langevin equations with a non-conserved 
order parameter, the equilibrium dynamic scaling fixed 
point remains stable, and governs the asymptotic critical 
behavior. Thus, the ensuing static critical exponents are 
those of the n-component Heisenberg model, accompanied 
with the equilibrium SSS strong-scaling dynamic expo- 
nent z = d/2 (at least to one-loop order). Novel dynamic 
and static scaling exponents are found in the extreme sit- 
uations t^/t'i = o or oo only, and may be related to 
the results for the isotropic non-equilibrium model stud- 
ied earlier ||. The latter are of course contained in our 
present more general study as a special case. These re- 
sults once more underline previous observations that non- 
equilibrium generalizations of dynamical models with a 
non-conserved order parameter asymptotically, i.e., in the 
vicinity of the critical point, display the scaling behav- 
ior of the corresponding equilibrium dynamic universality 
class. 

In Sec. ^, we consider the anisotropic non-equilibrium 
model J with conserved order parameter. In contrast to 
the anisotropic non-equilibrium SSS model, we shall find 
that only the spatial sector with lower noise temperature 
T 1 - < 7 1 " becomes soft at the transition, while the dii- 
dimensional sector remains uncritical. As a consequence, 
long-range elastic (uniaxial pseudo-dipolar for d± = 1) ef- 
fective interactions are generated, as in the two-tempera- 
ture model B . The equilibrium dynamic fixed point 
becomes unstable, and the phase transition is described by 
a novel universality class, characterized by reduced upper 
and lower critical dimensions for both statics and dynam- 
ics. For general dimensionality of the soft sector, we find 
runaway renormalization-group flow trajectories, and per- 
turbational methods appear to break down. Presumably, 
this indicates either strong-coupling scaling behavior, or 
perhaps even the absence of a non-equilibrium steady state 
at the critical point. However, we can obtain the new expo- 
nents in a double expansion about both the upper critical 
dimensions, and the dimensionality of the hard sector du 
(which amounts to an expansion about the equilibrium 
theory). Yet, as this expansion clearly becomes invalid at 
some critical value of djj < 1, at least to one-loop order, 
this formal expansion should not be taken too seriously 
for the description of a real physical system. At any rate, 
this model characterized by a conserved order parameter, 
and spatially anisotropic conserved noise, definitely leads 
to a novel dynamic universality class. Finally, in Sec. ^| 
we summarize our results again, and discuss their impli- 
cations. In the Appendix, we present the Ward identities 
stemming from the 0(n) symmetry of the non-equilibrium 
models under investigation here. 



2 Model equations 

In this section, we briefly outline the basic model equa- 
tions for our anisotropic, non-equilibrium generalization of 
both the SSS model Q and model J |ij . The equilibrium 
characteristics of these dynamic models were summarized 
at length in Ref. and we refer to this paper and the 
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original equilibrium literature (see Refs. |15|J16| , |17| , |18| 
for further details. We shall largely use the notations in- 
troduced in Ref. j5|, if not explicitly mentioned otherwise. 

We consider a second-order phase transition for an n- 
component vector order parameter Sq , a = 1, . . . , n (we 
denote unrenormalized quantities by the subscript "0"). 
As we furthermore assume isotropy in order parameter 
space, the static critical properties will be described by an 
0(n)-symmetric <fi Landau-Ginzburg- Wilson free energy 
in d space dimensions, 



H 



/( n 1 n 

d d x y$>o(x) 2 +2E[ V <?o(x 

n 



Uq 

"4! 



(2.1) 



where ro = (T — T®)/T® is the relative distance from the 
mean- field critical temperature T®. This effective free en- 
ergy determines the equilibrium probability distribution 
for the vector order parameter Sq , 



P cq [{sS}} 



-H[{SZ}]/k B T 



JV[{S«}]e- H ti s S}]/ k v T 



(2.2) 



Following standard procedures, one may then compute 
the two independent critical exponents, e.g., r\ and v, by 
means of perturbation theory with respect to the static 
non-linear coupling uq and by employing the renormal- 
ization group procedure, within a systematic expansion 
in terms of e = 4 — d about the static upper critical di- 
mension d c = 4. Here, rj describes the power-law decay 
of the order parameter correlation function at criticality, 
(S a (x)S l3 (x')) oc l/|x-x'| d - 2+? >, or, equivalently, of the 
static susceptibility, %(q) cx l/q 2 ~' n , and the exponent v 
characterizes the divergence of the correlation length as 
T c is approached, £ cx \T — T c \~ u . Notice that fluctuations 
also shift the true transition temperature T c downwards as 
compared to the mean-field critical temperature T°, i.e., 
r 0c = T c - T c ° < 0. 

In order to correctly describe the critical dynamics for 
an 0(n)-symmetric system, one needs to take into ac- 
count all the slow modes. Generally, in addition to the 
order parameter itself, these comprise the diffusive modes 
associated with the conservation law connected with the 
rotational symmetry. The SSS model thus consists of dy- 
namically coupled Langevin equations for a non- conserved 
n-component order parameter Sq and n(n — l)/2 con- 



-M, 



0a 



served generalized angular momenta M 
Physical realizations of this model are the critical dynam- 
ics of the XY model in — 2), also called model E 
with the order parameter components Sq and Sq, and 
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the conserved quantity Mg 2 = Sq, which generates rota- 
tions in the icy-plane; and the dynamic critical behavior 
of isotropic antiferromagnets (n — 3), known as model G 
with S'§, Sq, and Sq representing the components of 
the staggered magnetization, and Mg 12 = M§ , Mg 3 = Mg , 
and Mg 3 = — Mg denoting the components of the magne- 
tization itself, which are conserved and can be identified 
with the generators of the symmetry group 0(3). 

The variables Mq 13 represent non-critical quantities, 
and their coupling to the order parameter fluctuations Sq 
is of purely dynamical character. Hence it su ffices to sim- 
ply add a quadratic term to the hamiltonian ( |2.l| ), 

H[{S%}, {Mf}] = H[{S%}} + ljd d xJ2 M o Q/3 W 2 ■ 

a>/3 

(2.3) 

With this free energy funtional H , the coupled non-linear 
Langevin equations defining the SSS model read |l5| , [l6| ,p| 

dS$ _ v SH p 5H 



dt a "^8Mf 10 "SS"q 

= g J2 M o Ps o-Mr -V 2 )S% 



s U oa qP qP _1_ tot 



(2.5) 



and 



dM, 



a/3 



dt 



= ~9o 



SH p _ SH 



6H 

SMf 



n ot/3 



-g [SZV 2 S 3 -SSV'S l 



qa 



-D V 2 Mf 



n 



a0 



(2.6) 



(2.7) 



Here, go denotes the strength of the reversible, so-called 
mode-coupling terms, and £ a and rf 1 ^ represent fluctuat- 
ing forces with zero mean, (£ a (x, t)) = 0, (?7 a ^(x, £)) = 0. 
In order to fully characterize the dynamics, we further- 
more need to specify the correlations of these stochastic 
forces; for the order parameter, we simply assume a Gaus- 
sian distribution for the with the second moment 

(C (x, t)C V, 0) - 2A «5(x - x') 5(t - t') 6 a ? , (2.8) 

corresponding to non-conserved white noise. On the other 
hand, the conservation law for the generalized angular 
momenta Mq^, and the antisymmetry with respect to 
the tensor indices a, (3 implies that (r] a P(x, t)rf s {x!, t')) oc 
-V 2 <5(x - x.')5{t - t') (5 af) 5~< s - 6°* SPi). This functional 
form, prescribed by the Einstein relatio n en suring that at 
long ti mes t he equilibrium distribution ( |2.2[ ) with the free 
energy (2.3) will be attained, provides us with the possi- 
bility to impose a spatially anisotropic form of detailed- 
balance violation through the prescription 

(77^ (x, i W, 0) = -2 (SH V 2 + DjM) <5(x - x') x 

x6(t-t') (5 af3 5~< s -5 a5 5^) . (2.9) 



We may interpret this as follows. First, consider the 
special case D l l = = Dq, which is in fact the model in- 
vestigated in Ref. Q. We can then identify Ao/Ao = k^Tg 
and Dq/Dq = /cb^m as the temperatures of the heat baths 
coupling to the order parameter and conserved quanti- 
ties, respectively. For T$ = Tm, obviously detailed bal- 
ance holds. More generally, the ratio Tq = Tm/T$ de- 
scribes the extent to which this equilibrium condition is 
violated (in Ref. ||, 0q = 1/T was used instead); for 
Tq < 1, energy flows from the order parameter heat bath 
into the conserved-quantities heat bath, and vice versa 
for Tq > 1. Notice that as we are really interested in 
the vicinity of the critical point, Ts ~ T c , and Tq thus 
gives the heat bath temperature of the generalized angu- 
lar momenta measured in terms of T c . While new critical 
behavior ensues for either Tq = or Tq — 00, a one- 
loop renormalization-group analysis shows that for any 
< Tq < 00 the asymptotic critical properties are those 
of the equilibrium model, i.e., the renormalized T — ► 1 
under scale transformations In the anisotropic gen- 



eralization (2.9), we furthermore allow for different tem- 



peratures in du- and dj_-dimensional sectors in space (with 
d\\+d± = d), with associated distinct temperatures (again, 
essentially measured with respect to Ts ~ T c ), 



7t||/x _ 
1 o — 



ll/J 



Dq Ao 



(2.10) 



Therefore there appears an additional degree of freedom 
here, namely the ratio 



^0 = Tq /Tq 



(2.11) 



of the conserved-noise heat bath temperatures in the trans- 
verse and parallel sectors, respectively, which without loss 
of generality we may assume to be in the interval < Co < 
1, where erg — 1 obviously corresponds to the equilibrium 
situation. The one-loop renormalization-group flow equa- 
tions for this non-equilibrium SSS model with spatially 
anisotropic conserved noise will be derived and studied in 
Sec. |. 

Model J (in the nomenclature coined in Ref. O) cor- 
responds to the situation where the conserved quantities 
associated with the O(n) symmetry are in fact identical to 
the order parameter fluctuations themselves. The physi- 
cal realization, with n — 3, is the critical dynamics of 
isotropic ferromagnets, with Sq, Sq, and Sq denoting the 
three components of the conserved magnetization vector 
[ [l9|]2"o[| . The three coupled Langevin equations for model 
J read 



osi 

dt 



-^E^ 7 S^ + AoV 2 || + r (2.12) 
0,7 SSq ob 



-3o £ e^S^SZ + A V 2 (r - V 2 )S « 

/3,7 

^0 ^2 Q a 

6 



(2.13) 
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where go again denotes the reversible mode-coupling con- 
stant, and (C"( x i t)) — for the stochastic forces. As now 
the order parameter is conserved, we may introduce spa- 
tially anisotropic noise in the correlat or of the associated 
conserved noise, in analogy with Eq. ( |2.9| ) for the gener- 
alized angular momenta of the SSS model, 

(C a (x,t)C /3 (x',t')> = -2 (AJiVjj + A ^Vi) «5(x-x') x 

x5(t-t')5 af3 . (2.14) 

In the corresponding isotropic variant where Aq = Xq, a 
simple rescaling of the non-linear couplings uq and go can 
absorb the effects of detailed- balance violation entirely; 
this demonstrates immediately that this model is asymp- 
totically governed by the equilibrium critical exponents 
H. Here, on the other hand, we may define different ef- 
fective temperatures for the longitudinal and transverse 
sectors, respectively, 



in a probability distribution for the dynamic fields S§ and 



T, 



I /J 



III/- 1 1\ 



(2.15) 



and thus once again a novel variable emerges, namely the 
ratio cr = A^/AjJ = T^/T^j, with < a < 1. The dra- 
matic implications of this spatially anisotropic conserved 
noise will be investigated in Sec. ||. 



3 The anisotropic non-equilibrium SSS model 

We now proceed by considering the critical behavior of 
our anisotropic non-equilibrium version of the SSS mode l 
|PJ||, as defined by the Lan gevi n equa tion s (2.5) and (2.7), 



with the noise correlators (2^) and (2J3), respectively. In 



Sec. 3.1 , we perform the perturbational renormalization to 
one-loop order, following the procedures that were already 
employed in Ref. J5|. These in turn constitute the appro- 
priate generalization of the equilibrium renormalization 
scheme, see Refs. 



16,221]. From the renormalization con- 



stants (Z factors) that render the field theory finite in the 
ultraviolet (UV), one may then derive the RG flow func- 
tions which enter the Gell-Mann-Low equation. This par- 
tial differential equation describes how correlation func- 
tions change under scale transformations. In the vicinity 
of an RG fixed point, the theory becomes scale-invariant, 
and the information previously gained about the UV be- 
havior can thus be employed to access the physically in- 
teresting power laws governing the infrared (IR) regime at 
the critical point (t cx T — T c — > 0) for long wavelengths 
(wavevector q — > 0) and low frequencies (uj — > 0) . 



3.1 Perturbation theory and renormalization 

3.1.1 Dynamic field theory 

As a first step, we tr anslate the Langevin equations ( |2.5| ) 
and (2/?), with (2^8) and (2JS), into a dynamic field the- 
ory, following standard procedures 0,Eq,E|. This results 



Mf, 



P[{SZ,Mf}) cx / V[{iSg}] I V[{iM^}} 



x e 



J[{SS},{SS},{M^},{M^}] 



(3.1) 



with the statistical weight given by Janssen-Dc Dominicis 
functional J — Ji iar + Jiel + Jmc- Its harmonic part, in 
terms of the original dynamic fields So(x,t), Mq (x, t) 
and the auxiliary fields Sq (x, t) , (x, t) reads 

J haI [{S^},{S^},{Mf},{M^}} = 
= Jd d x Jdt^XoS^- 
d 



Of 



+ A (r - V 2 ) 



°0 



A I, 



a>/3 

e^(|-w 2 )m ^| 



(3.2) 



while the static non-linearity leads to a relaxation vertex 
■W^MS?}] = -Aoy J d d x J dt^SSS^S^Si; , 

a. (3 

(3.3) 

and the purely dynamic couplings generate the mode- 
coupling vertices 

Jmc[{Sn,{S^},{K },{M^}} = 

r d d x JdtY.lssK'sg- 



9a 



a,/3 



-iM ^(5 Q V 2 ^-^V 2 ^)| . (3.4) 

As usual, the harmonic part (|]^) defines the propagators 
of the field theory, while the perturbation exp ansi on is 
perf ormed in terms of the non-linear vertices ( |3.3|) and 
(3.4). Notice that the existence of the reversible forces 
( 3.4) does not show up in dynamic mean-field theory (van 



Hove theory), which in fi eld-theory language is based on 
the harmonic action (3.2) only. 

We can now construct the perturbation expansion for 
all possible correlation functions of the dynamic and aux- 
iliary fields, as well as for the associated vertex functions 
given by the one-particle irreducible Feynman diagrams. 
A straightforward scaling analysis yields that the upper 
critical dimension of this model is d c — 4 for both the re- 
laxational and the mode-coupling vertices. Therefore, for 
d < 4 the perturbation theory will be IR-singular, and 
non-trivial critical exponents will ensue, while for d > 4 
the perturbation theory contains UV divergences. In or- 
der to renormalize the field theory in the ultraviolet, it 
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suffices to render all the non- vanishing two-, three-, and 
four-point functions finite by introducing multiplicative 
renormalization constants. This is achieved by demanding 
the renormalized vertex functions, or appropriate momen- 
tum and frequency derivatives thereof, to be finite when 
the fluctuation integrals are taken at a conveniently chosen 
normalization point, well outside the IR regime. We shall 
employ the dimensional regularization scheme in order to 
compute the emerging momentum integrals, and choose 
the renormalized mass r = 1 as our normalization point, 
or, sufficient to one-loop order, tq = ro — r$ c — /j 2 . Notice 
that /i defines an intrinsic momentum scale of the renor- 
malized theory. The Gell- Mann-Low equation can subse- 
quently be used to explore the dependence of the renor- 
malized correlation or vertex functions on /z, and thereby 
obtain information on their scaling behavior. 



point, which yields after evaluating the integrals in dimen- 
sional renormalization |23j| , 



{Z~Z S 



1/2 = 1 



(3.10) 



d ll T ll 



(l + ^o) 2 



Here, e = 4 - d, A d = T(3 - d/2)/2 d - 1 7 r d / 2 is a geomet- 
ric factor (non-singular n ear d c = 4), T " and To 1 - denote 
the ratios defined in Eq. (2.10), and the effective dynamic 
couplings are 



Wq 



Ao_ 
Do 



fo = — fo = — ,9 ° 
Ao Aq Aq-Dq 



(3.11) 
(3.12) 



3.1.2 Dynamics: Vertex function renormalization 

The UV-divergent two-, three-, and four-point vertex func- 
tions or their derivatives that require multiplicative renor- 
malization are 9 w r Q ~ M (q, w), d q 2T Q ~ M (q, u), and d q 2 
r o5iEr( < J» w )i r oss( < J» w )> 9 u r ~ s (q,uj), ^/^(q.u;), 
and d q 2r ~~(q,uj); r Q ~ SM (-q, -u; § - p, f; § + p, f) 



and d( n .r,\r r 



IP) OMSS 



(-q, -w;f-p, 



2 ' 2 



p, and, at last, 



T ~ sss (-q, -u; §, f ; § , f ; f , f ). On the other hand, we 
have four fluctuating fieldsJM^, M Q/3 , S$, Sft) and the 
seven parameters Do, Dq, Ao, Ao, to, ffoj & n d uo available; 
this leaves us at liberty to choose one of the renormaliza- 
tion constants in a convenient manner. 

Starting with the two-point functions T ~ M (q, ui) and 
T Q g S (<li w) for the conserved quantities and order parame- 
ter fluctuations, respectively, we immediately note that as 
a consequence of the momentum dependence of the mode- 
coupling vertices 



d{iw) 



[q = 0,u) = 1 



(3.5) 



to all orders in perturbation theory. Upon defining renor- 
malized fields according to 



M 



«/3 - 7 1 ' 2 M a ? 



Ztl'M, 

M 
rl/2 oa 



> 



M a ? = Z\1 2 M? P 



s ' — s 



(3.7) 



and using T~ M = (Z~Z M ) 1/2 D njC?M , we thus obtain 



the exact relation 



Z~Z M = 1 



MM ' 



(3.8) 



At this point we utilize our freedom of choice to set 



Z 



M 



J M 



1 



(3.9) 



Similarly, for the order parameter fields we demand 
that dr~ s (q = 0,uj)/d(iu)) be finite at the normalization 



Notice that in equilibrium, when Tq = T^ = 1, Eq. ( (3 . 1 C ) 



yields Z~Z 



1/2 



S^S 



case, with T " 



= 1 to this order. In the isotropic noise 
Tq i , or equivalently, either dn or d± = 0, 



we recover the result cited in Ref. |5J] . In the same way, all 
the subsequently found Z factors reduce to the results for 
the isotropic non-equilibrium SSS model when Tq = Tq 1 - = 
To, and these in turn to the well-established equilibrium 
expressions for To = 1. 

As a next step, we compute the three-point vertex 
function f ~ M at zero external momenta and frequen- 
cies. Upon denning the dimensionless renormalized mode- 
coupling constant 

g = Z^g A 1 J 2 ^, 
this provides us with the product of Z factors 

(Z~Z S Z M Zg)^ 2 = 1 + 

wo/oAjM" 



(3.13) 



(3.14) 



d\\ „n d i „ i 

i - -At* - ^ 



(l + ^o) 2 

Direct comparison with Eq. ( 3.1 0| ) implies 



z g — z m — Z M ~ 1 > 



(3.15) 



where the choice ( |3.9| ) was employed. As shown in the Ap- 
pendix, as a consequence of the 0(n) invariance and the 
(3-6) fact that the Mg are the generators of the rotation sym- 



metry group, one may derive a Ward identity leading to 
the exact relation Z 3 Zm = 1- In equilibrium, this result is 
trivial, and a simple consequence of the fact that the con- 
served fields Mq 13 are non-critical. The absence of field 
renormalization therefore follows directly from the purely 
quadr atic appearance of the fields Mq in the hamilto- 



nian ( |2.3| ), which immediately implies that the static re- 
sponse function for the generalized angular momenta is 
Ao(q, u) = 0) = 1 exactly. While this relation holds even in 
our varian t of the SSS model with dynamic anisotropy, see 
Eq. ( 3.26| ) below, one cannot directly infer the renormal- 
ization constant Zm therefrom, as there is no fluctuation- 
dissipation theorem to relate this response function with 
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the corresponding correlation function in the general non- 
equilibrium situation. Therefore Z g needs to be computed 
explicitly from the three-point vertex function, or inferred 
from the above-mentioned Ward identity. 

The other vertex function renormalizing the mode- 
coupling constant g is r Q ~ ss (-q, -ui; f-p, ~ ; § +p, ~) 
= 2go(q-p) + 0(<7o)- Thus, in order to obtain Zs, we need 
to take a derivative with respect to (q • p) ; but owing to 
the dynamical anisotropy, the result depends on whether 
the components % and Pi lie in the longitudinal or trans- 
verse sector in momentum space, respectively. This means 
that we have to introduce different field renormalizations 
Zg and Zg for the longitudinal and transverse field fluc- 
tuations, for which one then finds to one-loop order 



Z. 



I/-L 



1 

n- 1 



(3.16) 



1 



2c 
n — 1 
~6e 



d ll T ll 









-t n 



foA d n- e 
(l + wo) 2 
2w Q 



(l + wo) : 



Obviously, for T Q = Tq 1 , this novel distinction between Z s 
and Zg disappears, and the isotropic result of Ref. ||] is 
recovered. Combining Eq. ( 3. IE ) with Eq. ( |3.10|) yields 



zi /J - = 1 
s 

n — 1 



2c 



1I T II 



4 -L r^J. 



±- 



1 dj 



(ic 



1 - Aw 
(l + wo) 2 
1 + 2w 



(3.17) 



(l + wo) 5 



foAd^" 



Next we define renormalized transport coefficients and 
noise strengths according to p4| 



DllA = ^|/- L £) , (3.18) 



A" 



A" 



7 II/-L x 



(3.19) 



where we allow for different renormalizations in the paral- 
lel and transverse sectors of the originally isotropic param- 
eters D , A , and A . The renormalized noise coefficients 
can be obtained by demanding that the vertex functions 
d i\ /± . r MM^WI^ u) = °)U|iM=o and I~~(q = 0,^ = 0) be 
UV-finite. This yields, with Z~ = 1, 



IIA 



and 



D 



n - 1 



M 

1 1_ foA d n~ 
+ 2e T \\/± 



o|| T || flx-i 



1 + U) 



(3.20) 



(3.21) 



While this product is still isotropic, the anisotropy in the 
field renormalization (3.17) induces different renormalized 
order parameter noise strengths in the longitudinal and 
transverse sectors. 



In the same manner, from d q %T~ M (q> w = 
d q ? r~ s (ciy ± ,u = 0)| q||/J = o we obtain 



Zd = 1 + Y e f^ A <i^ 



and 



(z~z s fi 2 z\^ = i + 



n - 1 f Q A d fi~ 



> S -oj ^\ 

n — 1 



1 - -4T " - -^T, 



1 + w 

wofoA d H~ 



< \ <T U d °J (l+w ) 3 
a - 1 d±/\\ f_w _ i \ wg/ yl d /i _£ 



1= o and 



(3.22) 



(3.23) 



=F 



3e 



(l + ^o) 3 



Thus, to one-loop order at least, we need not distinguish 
between the renormalized diffusion constants Z?" and D ± . 
This concludes our multiplicative vertex function renor- 
malization for the dynamical parameters of the non-equili- 
brium SSS model with dynamical anisotropy. Notice that 
the anisotropic parts of the renormalization constants are 



1 ■ 



For T, 



il _ 







rpA. _ 1 



always proportional to T 
we recover the equilibrium results where Z~ = Zjj and 
Z~ = Z\, reflecting the Einstein relation. 



3.1.3 Statics: Response function renormalization 

In order to define the "static" limit of the intrinsically dy- 
namic model under consideration here, we compute the re- 
sponse functions for the generalized angular momenta and 
the order parameter components, and then take the limit 
uj — > there. By adding external fields to the hamiltonian 



(2.3), one may show that the dynamic susceptibilities for 
the conserved quantities Mq and the order parameter 
fluctuations Sq are given by 

X (q,w) = J T ~ M (-q,-^)- 1 x 

x D q 2 + 2 9a r Q ~~ S] (-q, - w)] ,(3.24) 
Xo(q,w) = /^(-q,-^ 1 x 

x Ao-go^-j-^-q,-^) , (3.25) 

respectively pTp| . Notice that composite-operator vertex 
functions enter these expressions, which in general implies 
that new renormalization constants are required to remove 
the UV singularities of the response functions (equiva- 
lcntly, one may utilize the Z factors obtained from the 
multiplicative renormalization of the vertex functions plus 
appropriate additive renormalizations pl|,p[). 

Yet one may show to all orders in perturbation theory 
that 

F omaM^) = iuJ + D °<? 2 + 2 5o r f;|s S1 (q> w ) > (3-26) 
and consequently 

X o (q,o; = 0) = l , (3.27) 



8 



Uwe C. Tauber et al.: Non-equilibrium critical behavior of 0(n)-symmetric systems 



which means that there is no additional renormalization 
required here. On the other hand, the static limit of the 
order parameter susceptibility is in fact singular, which 
leads us to define the corresponding renormalized response 
function via 

X (q,o;) = ZXo(q,w) . (3.28) 

The new renormalization constant Z is determined by de- 
manding that £^2 x(q||/j_, w = 0) _1 | q||/i= o be UV-finite; 
once again, the result is different in the longitudinal and 
transverse momentum space sectors: 



= i 



71-1 



c 

Id 



1 - — r?n - 



*-L t-i-L 



w foA d p~ 



(l + ^o) 3 



=F 



-L/ll 



3c 



(l + ^o) 3 



(3.29) 



To one-loop order, Z^ 1 - as well as Z'j/ and Z~ , do not 

contain the static non- linear coupling uq- 

This is different for the remaining two Z factors needed 
for the renormalized dimensionless distance from the crit- 
ical point t and the static non-linearity u ]24| , 



r" 



= Zy ± T fl- 2 , 

= Z^ 1 - u A d [i 



to = r - r 0c , 



(3.30) 
(3.31) 



The fluctuation-induced T c shift is determined from the 
criticality condition xo(q = 0,ui = 0) _1 = at ro = rg c 
(tq = 0) with the result 



roc = 



■ u 



1 



k roc - 



+ (n- 1) X 



(3.32) 



fa 



1 



kWo(r 0c + k 2 ) + k 2 



In principle, these momentum integrals should be evalu- 
ated with a finite upper cutoff, which underlines the non- 
universality of the T c shift, i.e., its dependence on short- 
distance properties. However, if we choose to evaluate the 
momentum integrals by means of dimensional regulariza- 
tion, we are led to 



2A ri 



Foe 



(d-2)(4-d) 



■U 



n - 1 



w 



(3.33) 



w 



1 + Wq 



dn 



where wc have defined 



i-^rj-^],., 



A 



Uo = — u 



(3.34) 



Rendering x(q = 0, ui = Q) 1 UV-finite, after substituting 
To — tq + roc, then yields the isotropic product 



Z Z T = 1 7r — u A d fi 



Gc 



(3.35) 



71-1 



1 _ _i T ll _ ^ T< 



wofoA d n 

(i + w a y 



The anisotropy in the Z 11 ^ from Eq. Q3.29Q then induces 
the different Z factors for the longitudinal and trans verse 
momentum space sectors anticipated in Eq. ( 3.30| ). We 



note that alternatively we could have used 7—^,(0, 0), pro 



viding the (isotropic) combination (Z-Zs) 1 ^ 2 Z\Z T . (In- 
deed, the anisotropic contributions to Z^ 1 - [Eq. ( |3.29 )] 



and to [Eq. fl3.23|) 1 are identical.) 



Finally, we need the renormalization of the static cou- 
pling Mo, which we may obtain from the four-point func- 
tion r~ sss at vanishing external wavevectors and frequen- 
cies. To one-loop order, there appear ten Feynman dia- 
grams, the contributions of nine of which cancel in equi- 
librium (and only there!). A somewhat tedious calculation 
eventually gives 



(Z-Zs) 1 ^ 2 Zs Z\ Z u 
1 — u A d fi 



(3.36) 



1 IoA d [i~ 



Gc 



71 



1 



e 

n- 1 



II rp\\ 
~d ° 



C 



1 



J- n 



w 

foA d n~ 



1 _ *MI _ 



d " I {l+wo) 2 



(1 + w )u 



The anisotropies in zV^zV 1 - again cause the differences 

in Zu and Z^. No further UV renormalizations are re- 
quired, and we may now turn to the analysis of the ensuing 
RG flow equations. 



3.2 Discussion of the RG flow equations 

3.2.1 RG equations for the vertex and response functions 

By means of the above renormalization constants, we can 
now write down the RG (Gell-Mann-Low) equations for 
the vertex functions and the dynamic susceptibilities. The 
latter connect the asymptotic theory, where the IR singu- 
larities become manifest, with a region in parameter space 
where the loop integrals are finite and ordinary "naive" 
perturbation expansion is applicable, and follow from the 
simple observation that the "bare" vertex functions do not 
depend on the renormalization scale /x, 



dp 



r. 



0M k S r M'S' 



({q,o;};{ao}) = 



(3.37) 



with {a } = g ,Dl,D^, £> ,Ao,Ao,t ,wo- For the non-equi- 
librium SSS model with dynamical anisotropy, we have to 
treat the longitudinal and transverse sectors in momen- 
tum space separately, i.e., we need to understand the flow 
of the renormalized set of parameters {a"} = g,D»,D», 



All,All,- 



and {a^} = g,D^,D^, X ± ,X ± ,t j 



re- 



spectively. Replacing the bare parameters and fields in 
Eq. ( |3.37 ) with the renormalized ones, we thus find the 
following partial differential equations for the renormal- 
ized vertex functions in the longitudinal and transverse 
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sectors, 
d 



{all/- 1 -} 



+ 2 Q s + 2 Cs 



Upon introducing renormalized anisotr opic co unter 
parts for the effective dynamic couplings ( |3.1l| ), ( |3.12 ) 
and (pTTo| ), 



xr M^M< Ss (M^IIA^M* 11 ' }J =0 • (3.38) 
Here, we have introduced Wilson's flow functions 



and 



C ll/X - a— 



AM 1 - = u — 



II/- 



lnZ 



S ' 



In 

o a o 



(3.39) 
(3.40) 

(3.41) 



w \\/± = 
/1IM = 

T ll/-L = 



AHA 

DllA 
AII/-L 



AII/-L AllADHA 
SUA A"A 



(3.47) 
(3.48) 
(3.49) 



£>"A AHA ' 

the zeta functions to one-loop order read explicitly 



H A _n-lf dj| || _ dx„ A (1-WIA)/IIM 
; s ~ 2 V d d / (1 + u.llA) 2 



(the index "0" indicates that the renormalized fields and 
parameters are to be expressed in terms of their bare coun- 
terparts prior to performing the derivatives with respect 
to the momentum scale //). Note th at C,~ = Cm = and 
( g = — e/2 as a consequence of Eqs._( 3.15 ) and ( 3.13 ). 



The Gell-Mann-Low equation ( 3.38 ) is readily solved 
with the method of characteristics fi — > fx£; this defines 
running couplings as the solutions to the first-order differ- 
ential RG flow equations 

^^^=CF ± Wa ll/± W, allA(l)=a«A. (3.42) 



The solutions of the partial differential equations ( 3.38 ) 
then read 



r 



M k S r M'S l 



(M,{qn/±,^};{a |IA } 



(3.43) 



exp 



rd /X (£') + 3^(0 



d(' 



xf 



M k S r M l S' 



In the same manner, one can solve the RG equations for 
the dynamic susceptibilities, with the results 



(3.44) 



X(M,{q|| /x ,u,};{a IIM }J = 
= x(jd, {q|| A ,a,};{all/ x W} 



and 



x (^ ; {qn/i^};{a" A } 

xx^Mw/xM-M'Hf)} 

where, in analogy with Eq. ( |3.40 ), 

C IIA =/x#- lnZNA. 



(3.45) 



(3.46) 



1 d i 



± 



6 


d I 




n — 1 j 




2 l 


n — 


l di/n . 


6 


d 1 




-(n-l) 


n — 


1 dx/|| , 



]_ l[yl ij 1 

d d 



(i + 2mJI/- l )/H/- l 
(l + iullA) 2 

/ll/J 



(l + lollA)2 

(l + 2iull/- L )/ll A 



(3.50) 



± 



d|| T || _ d± T 
d 



( 1+W ||M)2 

w \\/Xf\\/± 



(3.51) 



yll _ j 1 - 1 



d " y (i + W n/-L) 3 



=F 



(l + TfillA) 3 



(3.52) 
(3.53) 
(3.54) 



3 d 

C =-- 
9 ~ 2 ' 

d 2 Tll/ J 

n-l / _ dj| „ _ dx A (1 + 6 W HA)/1IA _ 
2 V d d / (1 + wHA) 2 + 



=F- 



n 



1 d. 



/II 



ah/ 

/•II/ 
V\ 



6 d 
1 

~ ~2 



(l + 2wllA)/HA 

(1 +W IIM)2 



-(n-l) 



,f" 



l + tollA 



(3.55) 
(3.56) 

T 



=F- 



(l + u-HA) 3 

-l^/ll^ii T x U^ ll/± ) 2 / ll/± (3 57) 

3 d V 1 J (l + u;IIA)3 ' (d - 57j 



CjA = -2- 



n + 2 
6 



5«A_ 



.(„-i) (i-^LTH-%r- , , 

• d d J (1 + iullA)3 



,IIA 2 /H 



_I _^/l f T ll _ t 



,11/- 



7" 



(1 + U-H/-L) 3 



(3.58) 
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C IIA = _ £ 



n- 1 



IIA 



8 



(3.59) 



d d 



±- 



-(n-l)(l- 

- 1 



6 



d 



I rj^\\ d_\_ rjnA 



(3 + 5^11/^)/"/^ 
(1 + ^IIA) 3 

6(7"^) 2 
1 + toll/- 1 ) ull/- 1 

(1 + 3^11/^ +4wll/ ±2 )/1l/ J 
(1 + wllA) 3 



± 



These results enable us now to study the scaling behavior 
of the non-equilibrium SSS model with dynamical noise 
in the vicinity of the different RG fixed points, which are 
given by the zeros of the appropriate RG beta functions 



(3 V = fi 



d_ 

Of! 



According to 



(3.60) 
(3.61) 



these govern the flow of the effective couplings T"'- 1 -, 
/I'/ -1 , and u"/- 1 , etc. under scale transformations [i — > fi£, 
and the fixed points {v*} where all f3 v ({v*}) = thus 
describe scale-invariant regimes. 



3.2.2 RG fixed points and critical exponents 

We begin with considering the RG flow of the anisotropy 
parameter 

a = T^/T 11 , (3.62) 

denoting t he ra tio of the effective conserved noise tem- 
peratures ( |3.49] ) in the transverse and longitudinal sectors 
(0 < a < 1). Obviously, a = 1 describes the isotropic 
fixed point. In order to assess its stability against the 
anisotropic non-equilibrium perturbation, we consider the 
RG beta function 



d_ 



(3.63) 



As to first order w" = w 1 - — w and f ± — — f (in 
the vicinity of the isotropic fixed point, this holds even 
beyond the one- loop approximation), we may write 



0a 



-<7(1-CT)/T« 



1 



1 1 + 3w 



2T\\T X 



6 (1 + w) 3 



(3.64) 



The expression in square brackets is positive, and thus, as 
to be expected on physical grounds, there are only two 
fixed points, namely a* = 1 and er* — 0, realized for 
T 1 - = and T" = oo. (Of course, if we allow for a > 1 
as well, there is also the fixed point a* = oo, realized for 
T" = and T x = oo; yet clearly the regimes a < 1 and 
a > 1 map onto each other through simply relabeling the 




Fig. 3.1. Equilibrium (center) and non-equilibrium (corners) 
fixed points of the non-equilibrium SSS model with dynamical 
anisotropy in the [T [[ /(l + T 11 ), /(I + T x )]-plane. 



indices ||<->_L.) Furthermore, in the 1R regime (£ — > 0), 
if < er < 1 initially, (3 a < 0, and thus a(£) grows un- 
til it reaches the isotropic fixed point a* = 1 (and con- 
versely, if 1 < er < oo, then (3 a > and a(t) decreases 
towards a* = 1). Thus, the isotropic fixed point is stable 
against the spatially anisotropic perturbations in the noise 
correla tor of the conserved generalized angular momenta. 
Figure 3.1 depicts these various fixed points in the non- 
equilibrium SSS model with dynamical anisotropy, and the 
parameter flow in the [TH/(1 + TH), T ± /{1 + T^-plane. 
The center of this diagram represents the equilibrium SSS 
dynamic fixed point (cr = 1, T = 1), and is attractive for 
the RG flows originating from any point inside the de- 
picted square. Precisely on the edges in parameter space, 
we find the two isotropic non-equilibrium fixed points with 
<T = 1 and T — in the lower left, and T = oo in the upper 
right corners, respectively, and the two anisotropic fixed 
points with a = and a = oo in the lower right and upper 
left corners. Notice that the RG flows on this critical sur- 
face in parameter space which start in the vicinity of these 
anisotropic non-equilibrium fixed points tend towards the 
isotropic non-equilibrium fixed points, but eventually end 
up at the isotropic equilibrium fixed point provided that 
initially < TH < oo and < T 1 - < oo. 

Before we investigate the properties of the anisotropic 
fixed point a* = (T 1 - = 0, = oo), let us briefly sum- 
marize the behavior of the isotropic model with a* = 1. 
(For more details, and for a graph depicting the vari- 
ous equilibrium and non-equilibrium fixed points in the 
isotropic parameter subspace, we refer to Ref. ]5[].) Set- 
ting TH = T 1 - = T in the flow functions $.hty-<$.S§), we 
find 



f3 T =T{(~-( D +C 



Cl = 



(3.65) 



1 , .l + 7w + Aw 2 

h [fl — 1) ; 

T ' (1 + u;) 3 



Clearly, the only possible fixed points here are T* = 1, 
T* = 0, and T* — oo. At the equilibrium fixed point 
T* = 1, one finds Ca = ~{n- l)//(l + iu) and ( D =-f/2, 
and the beta functions for the couplings w — X/D and 
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/ = g 2 /XD read pUT 



Pw = w (Ca - Cd) 



-wf 



n-1 1 
1 + w ~ 2 



P f = / (2C 9 - Ca - Cd) = 

'n-1 1 
1 + w + 2 



= / 



-e+/ 



(3.66) 



(3.67) 



The IR-stable fixed point (to one-loop order at least) turns 
out to be the strong-scaling SSS fixed point 

w*=2n-3, f* = e + 0(e 2 ) (3.68) 

with equal time scales governing the critical slowing down 



of the order parameter (t$ oc £ 2 



,v ) and conserved 



generalized angular momenta fluctuations (t c M oc £ 



~ ZMV )^ respectively, 

z s = 2 + Cx , z m = 2 + Cd 



(3.69) 



Inserting the fixed-point values (3.68) yields the dynamic 
critical exponent 



d 



ZS = Z M 



(3.70) 



which is actually an exact result, pro vided z$ — zm and 
< /* < oo is finite, because Eq. ( 3.67| ) then requires 
that 2z = 4 + Q + Cd = 4 + 2 C g = 4 - e = d |7||l6). The 
equilibrium static critical behavior is described by the zero 
of the beta function 



n + 8 



(3.71) 



which yields of course the 0(n) Heisenberg fixed point 



6 



u 



e + 0(e 2 ) 



(3.72) 



n + 8 

with the associated two independent critical exponents 

n + 2 



v 



T) 



-1 



-c; 



n + 8 
-C* =0 + O(e 2 ) . 



(3.73) 
(3.74) 



For the isotropic non-equilibrium fixed point with tem- 
perature ratio T* — oo, the appropriate effective mode- 
coupling constant becomes 



j = T J= D D_l_ 
1 w J X A XD 



in terms of which the beta functions read Q 



Pw =w(C\- Cd) = -{n - 1) 



/^ = /(2C 3 + C 5 -2Ca-Cd) = 



= /(^-e + 2(n-l) 

Pz = u(cu + c~ x -a) 

n + 8 



w 3 f 
(1+w) 3 



(3.75) 
(3.76) 

(3.77) 
(3.78) 



~ 0/ n w(l + 3w + w 2 ) - 

-e + — — u - 2 n - 1) / 

(1 + w) 6 



Thus, the RG fixed points to 0(e) governing this scaling 
regime are 
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2(n- 1) 



(3.79) 



Notice that the fixed point o f the static coupling takes 
on twice the Heisenberg value (3.72); correspondingly, the 
"static" exponents will be modified as compared to the 
equilibrium results. E.g., the correlation length exponent 
now becomes 



i n + 2 e 

v =2-2 tC-tt 

n + 8 2 



0(e 2 



(3.80) 



instead of Eq. (3.74), while to one- loop order both the or- 
der parameter response and correlation function are char- 
acterized by the Wilson-Fisher exponent n = + 0(e 2 ) 
as in equilibrium ||. The characteristic time scales for 
the order parameter and the conserved quantities are now 
governed by different power laws, namely 



zs 



zm = 2 



(3.81) 



As T* — co means that effectively the heat bath for the 
conserved generalized angular momenta is at infinite tem- 
perature, there is effectively an energy current into the or- 
der parameter heat bath, but no feedback. This explains 
why we find the coupled SSS model dynamic exponent 
for the order parameter fluctuations, while the general- 
ized angular momenta correlations decay faster with the 
purely diffusive exponent zm — 2. Finally, at the critical 
point there are non-trivial noise correlations oc q p only for 
the order parameter noise, 



Ps = 2 + 



0(e 2 



while for the generalized angular momenta 



PM = 2 + C 



D 



(3.82) 



(3.83) 



as to be expected at infinite temperature Q. 

For the other isotropic non-equilibrium fixed point, be- 
ing characterized by T* — 0, the correct effective mode- 
coupling constant reads 



/' 



T \2 D XD 



(3.84) 



(called / in Ref. ||), and consequently (3 W = 0, which 
leaves the fixed point w* undetermined. The remaining 
RG beta functions are |fj| 



H =/'( 2 C 9 + 2Ct-3Ca-Ch) =f'[-£ + lf 



'/' J V a n s a ~™ ->dj - y ' ■ 2 

(3.85) 

= ! . (3-86) 
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with the O(e) fixed points 

r=2e + 0(e 2 ), u* = ^e + 0(e 2 ) . (3.87) 



As u* is identical to the Heisenberg fixed point ( 3.72[ ), it 
turns out that the static ex ponen ts are inde ed those of the 
equilibrium static theory, ( |3.73 ) and ( 3.74 ), both for the 



order parameter response and correlation functions. Now 
the energy current flows from the order parameter heat 
bath towards the conserved generalized angular momenta, 
and thus the mode-coupling effects become negligible for 
the order parameter fluctuations. The dynamic exponents 
therefore become model-A like, with 

z s = z M = 2 (3.88) 
to one-loop order, and with the critical noise exponents 



p s = 2 , p M = 2-e = d-2 



(3.89) 



Again, we observe that these anomalous power laws for 
the noise correlators apply for those quantities which are 
governed by the heat bath at lower temperature; here, the 
generalized angular momenta at effectively T = jf| . 

The above discussion of the isotropic fixed points fa- 
cilitates the interpretation of the novel fixed poi nts of the 
SSS model with dynamic anisotropy. From Eq. ( 3.64 ) we 
had already inferred that at the fixed point with a* = 0, 
one must have T" — oo and T 1 - = 0. In analogy with 



Eqs. ( 3.75 ) and ( 3.84 ), it is therefore convenient to intro- 
duce new effective mode-coupling constants in the longi- 
tudinal and transverse sectors, respectively, 



7" = ^ 7" 

.11 

Ti± _ f_ 

J rpj_ 



DW ff 2 



AH AIIDH ' 

„2 



2 

X 1 - D 



In terms of these couplings, one finds for a* 



T 



1 



Til 



(3.90) 
(3.91) 

(3.92) 



/ .11(1 + 7.11 +4.ll 2 ) + 

; V d ) (1 + .H) 3 

, n - 1 d± .11(1 + 3. 11 ) 



fox = -T 1 



3 d 
1 -T 1 - 



(l + .H) 3 



(3.93) 



+{ n-l)T- fl-^Tll V + J; X+ x t X) 
v ' \ d J (l + w-L) 3 

, n ~ l d \\ T \\ T ± 1 + 3. x ) 
3d (1 + .-L) 3 

which indeed lead to the expected fixed points 

T"* = oo , T^* — , (3.94) 



with Til* t 1 -* = o. 

In the longitudinal sector, we may then write the beta 
functions for w", /", and u" as 



fiwW — — 



n-1 



||4 7 



1 -I- 2— 

d J (l + .H) 3 



||3 7 || 

.11 / 



2(n- 1) / d|| 

3 l 1 + 2 T; (l + .ii) 3 



(3.95) 



(3.96) 



G 



d,l .ll(l + 3.H + .H ) 7 || 
1 ' d (l + .H) 3 ; + 



2{n-l)d± .II 3 / 11 



3 d (l + .H) 3 
yielding the stable fixed point 



II 

." = oo 



and consequently 

II * 

/ = 



2(n-l)(l+2d||/d) 
II „ _ 36dj|/d 



e + 0( e 2 ) 
+ 0(e 2 ) . 



(3.97) 

(3.98) 

(3.99) 
(3.100) 



(n + 8)(l + 2dy/d) 

For d±_ = 0, i.e., du — d, these expressions reduce to the 
isotropic fixed point (3.7E) with T* = oo, and in fact, we 
arrive at very similar results for the dynamic exponents, 



2 + C 



All 



= 2 + Ci» = 2 , 



(3.101) 
(3.102) 



as well as for the anomalous noise exponents, 



II „ 1 + 8d\\/d , 9n , 



2(l + 2d||/d) 



pL = 2 + c~ 



(3.104) 



compare Eqs. fDjl] ) and ( ||| ), ( |3~83| ). Moreover, we again 
find a non-standard correlation length exponent 



"Cn 



6(n + 2)d\\/d 
(n + 8)(l + 2d||/d) 



l-4d||/d 
"2(1 + 2d||/rf) 



£ + 0(e 2 ) 



(3.105) 



but, in addition, non-trivial Wilson-Fisher exponents de- 
scribing the critical decay of the order parameter response 
and correlation functions, respectively, 



T,ll = -C" = 



I* - 



1 - d\\/d 



'Is ~ — 



2(l + 2d||/d) 
1 - d||/d 



e + 0(e 2 ), (3.106) 



2(l + 2d|i/d) 



e + 0(e 2 ) ; (3.107) 
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remarkably, rji = —77" to 0(e). These anomalous critical 
exponents at the anisotropic fixed point, appearing in the 
longitudinal sector with effectively infinite heat bath tem- 
perature for the generalized angular momenta, are obvi- 
ously a consequence of the spatially extremely anisotropic 
noise correlations in the conserved quantities. They may 
perhaps be interpreted as remnants of the elastic, pseudo- 
dipolar interactions generated in a model with conserved 



order parameter and dynamical anisotropy |5||7 ,|8j • 

We now turn to the transverse sector, with the con- 
served quantities being effectively at zero temperature as 
T = 0. Quite as at the isotropic non-equilibrium SSS 
model fixed point with T* = 0, we find that (3 W ± = and 
hence w has no fixed value, while 



-e+——u 1 
6 



(3.108) 
(3.109) 



which yield the one-loop fixed points 

6 



f ,A - = 2e + 0(e z ), u A 



e + 0(e 2 ) . (3.110) 



4 The anisotropic non-equilibrium model J 

In this section, we study the critical behavior of our non- 
equilibrium version for model J (describing the dynamics 
of isotropic fcr roma gnets ) with dynamical noise, as defined 
through Eqs. ( 2. IS ) and ( 2.14 ). We start by computing the 
T c shift from the static susceptibility. As a consequence of 
the spatially anisotropic conserved noise with < , it 
turns out that the transverse momentum space sector with 
lower noise temperature softens first. Thus, at the critical 
point, the longitudinal sector remains uncritical ("stiff"), 
similar to equilibrium anisotropic elastic phase transitions 
[p5[ . It is then instructive to switch off the mode-coupling 
constant go, and first recapitulate the properties of the en- 
suin g tw o-temperature non-equilibrium model B ^|J7|^] . In 
Sec. O, we turn to the perturbational renormalization of 
the two-temperature non-equilibrium model J to one-loop 
order, and finally discuss the resulting RG flow equations. 



4.1 Dynamic field theory and the anisotropic T c shift 

The probability distribution for the dynamic field s Sf? 
(a = 1,2,3), equivalent to the L angey in equation (2. IS) 
with anisotropic noise correlator ( 2.14 ), reads 



The ensuing critical exponents are precisely those of the 
T* = isotropic non-equilibrium fixed point: 



P[{SS}}<* /P[{!5 a }]e J|{s »)' {s » }l , (4.1) 



and 







C T x. = 2 - 


n + 2 , , % 
, + 8 e + °^< 


(3.111) 


n L = 


-c J 


* X 

~ = Vs = - 


-cr = o + o( e 2 ) 


(3.112) 


4 


= 2 


+" Ca- 1 - = Z M 


= 2 + C D ± = 2 , 


(3.113) 


pi 


= 2 






(3.114) 


Pm 


= 2 




- e = d- 2 . 


(3.115) 



with the Janssen-De Dominicis functional J = Jhar+^rel" 
J mC) with the harmonic part 



J har [{5 «},{^}] = 



(4.2) 



We may interpret our results for the novel anisotropic the non- linear relaxation vertex from model B, 



fixed point as follows. For T = and T" = 00, the 
system breaks up into essentially independent sheets of 
dimension du with infinite heat bath temperature. The 
associated critical exponents are closely related to the 
isotropic ones at the fixed point with T* = 00. However, 
the additional d± dimensions a re refl ected in the anoma- 
lous Wilson-Fisher exponents ( |3 . 1 06 ) and (3.107), which 
are proportional to d± e (while the equilibrium exponent 
77 oc e 2 ). Fluctuations in the d± transverse direction are ef- 
fectively at zero temperature for the conserved noise, and 
consequently are governed by the critical exponents of the 
isotropic non-equilibrium SSS model with T* — 0. In the 
converse situation, with T" = and T 1 - = 00, clearly one 
must simply exchange the roles of the transverse and longi- 
tudinal sectors. Yet, we emphasize again that these novel 
fixed points with their associated rather bizarre critical 
behavior are unstable, and for any < 00 

initially, the static and dynamic critical properties of the 
system are asymptotically described by the equilibrium 
scaling exponents. 



J rcl [{5 "},{^}] = A ^ J d d x J dtJ^S^S^S^ , 

a, (3 

and the mode-coupling vertices from the spin precession 
forces, 

J mc [{5 Q }, {S%}} = ~g Q J d d x JdtY, ^SgSgV 2 ^ . 

j-a (44) 

In analogy with Eq. Q3.25j ) for the SSS model, the dy- 
namic susceptibility can be expressed as 

Xo(q,w) = ^-^(-q,-^)" 1 x 

x X q 2 +50-^5^] (-q,-^) • (4.5) 

From the ensuing expression (to one- loop order), we may 
determine the fluctuation-induced shift of the critical tem- 
perature. Because of the dynamic anisotropy appearing in 
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the noise correlator (2.14), however, the result depends on 
how the limit q ► is taken; upon defining q« — qcosO 
and q± — qsm.0, we find 



r Qc (0) 



a \\ T \\ 

i 



■ T J 



(4.6) 



9l 



6 U ° J k r 0c + k 2 2(d + 2)X 2 J k k 2 (r 0c + k 2 ) 

1 



1 



- t " cos^e 



T ± sin^ 



fi 



2{d- 



2) Ag J k k 2 (r 0c 



in contrast with the isotropic Eq. ( |3.32 ) for the SSS model. 
As T c = T c ° + roc, the phase transition will occur at the 
maximum of the function ro c ((9), which for T^ < Tq' oc- 
curs at — ir/2. The <ij_-dimensional transverse sector in 
momentum space thus softens first, and the true T c shift 
is given by 



roc = 



Til 



d \\ T \\ 



9l 



6 U ° J k roc + k 2 



2(d+2)X 2 J k k 2 (r 0c + k 2 ) 



, (4-7) 



or, after evaluating the integrals in Eq. ( |4.7| ) by means of 
dimensional regularization, as solution of 



2A rl 



Foe 



(d-2)(4-d) 



I rp\\ 

7 J ° 



d ±- rpl 



<i 



2{d- 



2) A 2 



(4.8) 



(We remark again, though, that a more physical way to 
compute this quantity would be by means of cutoff regula- 
rization.) For Tq = Tq- = To, we recover the equilibrium 
result with a rescaled coupling TqUq = XqUq/Xq, as to be 
expected (see Ref. || an d Eq. ( 3.33 ) for the SSS model 
with n = 3). Notice, however, that dynamical anisotropy 
(Tq =/= Tq 1 ), combined with the reversible mode-coupling 
terms, has a very drastic effect here: It renders the sys- 
tem soft only in the momentum subspace with lower noise 
temperature. This effect has a simple physical interpreta- 
tion: The T c shift is due to thermal fluctuations, which are 

reduced in the transverse sector (Tq 1 < Tq), and therefore 
lead to a comparatively stronger downwards shift in the 
longitudinal sector. 

In order to characterize the critical properties of our 
model, we may neglect terms oc qf, in the stiff momentum 

space sector, because Tq = r$ — tq c {0 = 0) remains posi- 
tive at the phase transition where Tq = ro — ro c (0 = it/ 2) 
vanishes. In analogy with the situation at anisotropic elas- 
tic structural phase transitions , or with Lifshitz points 
in magnetic systems with competing interactions pq , we 
thus have to scale the soft and stiff wavevector compo- 
nents differently, [q±] = fi, whereas [qii] = [q±]' 2 = /i 2 . 

Consequently, while [Xq] = /i° if we choose [u>] = /i 4 , 



we find for the longitudinal scaling dimension [A ] = [i~ 2 
which implies that the longitudinal noise strength becomes 
irrelevant under scale transformations. Allowing for dis- 
tinct couplings in the different sectors, one finds in the 
same manner that the ratios [Aq/Aq] = [AqMq/A^Wq 1 ] = 

[So/ffcH = I 1 2 au have negative scaling dimension. Thus, 
for an investigation of the asympotic critical behavior, 
the longitudinal parameters may be neglected as com- 
pared to their transverse counterparts, and can all be 
set to zero in the effective dynamic functional J. Upon 
rescaling the fields according to Sq — > (Xq / Xq ) 1 / 2 Sg , 
S% - (AqVA 1 ) 1 / 2 ^, defining 



Co 



M 



9o 




(4.9) 



and omitting the labels "_L" again for Ao and ro, the en- 
suing effective Langevin equation of motion becomes 



dsi 

dt 



= A 



Ao-rViSo ^2 s o s o 



coVf +Vi(r -Vi)| 5* " 



(4.10) 



/3,7 



with the noise correlator 

(C tt (x,<)c /3 (x',t')> ^-aAoVi^x-xO^-O^ ■ 

(4.11) 

These equations define the two-temperature non- equilibri- 
um model J. 

It is interesting to note that the anisotropy of the T c 



shift in Eq. (4.6) only occurs in the contribution oc g^. In 
the non-equilibrium model B with dynamical anisotropy, 
the criticality condition for the response function remains 
isotropic, at least to one-loop order. Thus, if one does not 
assume different critical temperatures in the purely diffu- 
sive non-linear Langevin equation to begin with, these are 
not generated, and one is not led to the two-temperature 
model B, which we shall discuss below, as the correct ef- 
fective theory for the phase transition. Instead, the non- 
equilibrium perturbations appear to be irrelevant to this 
order, and the model is asymptotically described by the 
equilibrium critical exponents of model B, i.e., the static 
exponents of the 0(n) Heisenberg model, accompanied 
with the dynamic exponent z = 4 — 77. 



4.2 The two-temperature non-equilibrium model B 

Before we turn to the analysis of the two-temperature 
model J, derived above as the effective critical theory for 
the non-equilibrium model J with dynamical anisotropy, 
we briefly summarize the results for the corresponding 
two-temperature n on-eq uilibr ium model B which 
is defined by Eqs. (4.1C) and (4.11) with vanishing mode- 
coupling term "qq = 0. We may thus generalize to arbitrary 
number of compenents n again. 
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With this simplification, the resulting purely relax- 
ational Langevin equation of motion can be written in 
the form 



dt 



. y2 6H eS [{S$}} 

AoV_l 



+ C 



(4.12) 



accompanied with the Gaussian noise ( 4.11 ). Notice that 
after the above rescaling, the Einstein relation between 
the diffusion constant and the noise strength is fulfilled; 
hence the two-temperature model B is effectively an equi- 
librium system, and describes diffusive relaxation into the 
stationary state with probability distribution P cq [{S , Q }] oc 
exp(— H c ff [{5q }]). The effective free energy here, 



1 f d d q v-^ coqj +q 2 j_(ro + q\) 



2 J (2tt) 



d ' 



^ Q (q)^o(-q) + 



' a,/3 



B> / M 
4! 



(4.13) 



contains long-range elastic interactions (uniaxial pseudo- 
dipolar for d ± — 1 |^7|), as is evident from the harmonic 
part of Eq. ( 4.13] ). These long-range, non-analytic inter- 
actions are generated by the dynamical anisotropy in the 
original non-equilibrium model ■ As critical fluctua- 

tions are now restricted to the e?x-dimcnsional subsector, 
one expects that the upper critical dimension of this model 
is reduced to 

df=A-d\\, (4.14) 

which is con firmed through a direct scaling analysis of the 
free energy ( [4.13 ). E.g., for a one-dimensional hard sec- 
tor (uniaxial system), the critical dimension above which 
mean-field exponents become exact, is d^} = 3. 

In order to compute the scaling exponents below the 
upper critical dimension, we have to renormalize the the- 
ory in the UV. To this end, we introduce renormalized 
fields and parameters according to 



S = Zl /2 S , 



A 



s 

Z\ A 



T = Z T T /X 



C : 
U 



Z c c , 
Z~u A(dn 



(4.15) 
(4.16) 
IM" £ , (4.17) 



where tq = tq — r$ c as usual, e = df — d = 4 — d — d\\ — 
4 — 2dn — d ± den otes the deviation from the upper critical 
dimension (4.14), and we define the anisotropic geometric 
factor 



A(d h d±) = 



r(3-d/2-d||/2)r(d/2) 

2<*-i 7r d/2 r(d±/2) 



A\/ 2 



(4.18) 



with A(0, d) = A d . 

Yet, these renormalization constants are not entirely 
independent of each other. First, as the two-temperature 
model B in the critical region is equivalent to an equi- 
librium system, there exists a fluctuation-dissipation the- 
orem that connects the imaginary part of the dynamic 



susceptibility xo (q, with the Fourier transform of the 
order parameter correlation function Co(x, i; x', t') 5 a/3 — 

(SS(x,t)Sg(x?,lf)), 



C (q,w) = - Imxo(q,w) 



(4.19) 



In terms of the two-point vertex functions, Cp( q, m ) = 

^^oss^'^^/l^oss^'^l 2 ' while for 9o = E( l- (0) re " 
duces to Xo(q.j w ) = ^oq±/r Q g S (— q, — w); thus the fluctu- 
ation-dissipation theorem can equivalently be written as 



2X q 2 j_ 



Imr ~ s (q,^) 



(4.20) 



Precisely the same relation must hold for the correspond- 
ing renormalized vertex function, which implies the iden- 
tity 

Z x = [Z s /Z~) . (4.21) 

Second, the equation of motion ( |4. 12| ) implies that the 
non-linear relaxation vertices are proportional to the ex- 
ternal momentum q\ , and hence the loop contributions to 
r" ~ (q, w) must vanish in the limit — > 0. Thus, to all 
orders in perturbation theory, 



r Q g S {q.\\,q± = Q,uj) = iu + \acoq\\ . 
This leads to the additional set of identities 



ZgZ s = 1 



Z\ z a 



1 



and thus, using Eq. (4.21), 

Zx = Z- 1 



Zs • 



(4.22) 



(4.23) 



(4.24) 



At last, because of the absence of the composite-operator 
vertex function in the relation between the dynamic sus- 
ceptibility and the two-point vertex function, we have 



with 



X(q,w) = Zxo(q,w) 



Z=Z S 



(4.25) 



(4.26) 



As in Sec. 3.2 for the non-equilibrium SSS model, we 
can now define Wilson's zeta functions via logarithmic 
derivatives of the Z factors with respect to the renormal- 
ization scale /1, 



Cs 

Ca 



d_ 
d_ 

d_ 

dfi 



lnZ~ 



IylZ s 



a o 



(4.27) 
(4.28) 
(4.29) 



where {a} = A, c, r, u, and write down the RG equations 
for the vertex and response functions, 

r^si (Mqih<uM;{a}) = 
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df 



xr~ ss (fd,{d h cu^h{a(e)}) , (4.30) 

x(Mq||,q-L,w};{a}) =cxp | - / ^^')^-| x 

xx(M^{q||,q±,w};{aW}) ■ (4.31) 

The general scaling form for the renormalized order pa- 
rameter response and correlation function thus obtained 
at an IR-stable fixed point u* reads pf 



/ \ -2-H7 ~ T 1\\ U 

\1± q ±- q ±. 



(4.32) 



C(r,c Lh q ± ,u>) = q ± >C I 

V 9x 9 -L ^-L / 

(4.33) 

where in addition to the usual static exponents r\ and v, 
and the dynamic exponent z, we have introduced a scaling 
exponent A originating in the intrinsic anisotropy of the 
system. Alternatively, we could have defined a set differ- 
ent transverse and longitudinal critical exponents, v±_ = v, 
z± = z, v\\ = u/(l+A), z\\ = z/(l +A) 0. Notice that as a 
consequence of the identity ( 4.26| ) , the exponent r\ governs 
the critical decay of the response as well as the correlation 
functi on, as is required by the fluctuation-dissipation the- 
orem ( 4.19 ). The critical exponents are readily identified 
with the fixed-point values of the flow functions 



n 



-Q, 
4 + Ca , 



A = 1 



-c 

c 

2 



(4.34) 
(4.35) 



The exact relation 
sequently 



con- 



Z\ = Z c 1 implies (\ = — Cc and 

l+A = z/2, (4.36) 

or 211 = 2, which reflects the mean-field character of the 
fluctuati ons in the stiff sector, and obviously holds when- 
ever Eq. ( 4.22 ) is valid. From the second equation in ( 4.24 ) 
we furthermore infer £x = Cs for the two-temperature 
model B, and hence M 



z = A-r) , A=l-r)/2 



(4.37) 



The task is thus to compute the remaining indepen- 
dent static exponents rj and v. To one-loop order, from 

h( d g 2 J 2 r ss( q W = 0c U w = )U-l=o one nnds z x = 
Zs = 1. Next, the criticality condition for the static sus- 
ceptibility yields the T c shift 



Foe 



u A(d\\,dx) 



3 (d + d ll -2)(4-d-d\\ 

(4.38) 

the denominator here indicates that in addition to the 
reduction of the upper critical dimension d^, the lower 



critical dimension appears to be lowered by d\\ as well 
to di c = 2 — du. Thus, in two dimensions, an ordered 
phase with long-range order may exist at low tempera- 
tures and is indeed found notice that the Mermin- 
Wagner-Hohenberg theorem is invalidated by the exis- 
tence of long-range elastic or pseudo-dipolar interactions 
in the system. Upon then replacing ro = To+ro c , rendering 
0, q_± 7 uj = Q)\q ±= o UV-finite gives 



Z T = 1 - 



(if 



■uoA(d\\,d±) fi~ 



(4.39) 
provides us 
(4.40) 



while the four-point vertex function F ' ~ . 

U o o o o 

with 

„ n + 8 „ . , , , . 
Z~ = 1 — u A(d\\ , d±_) (i 

Notice that the combinatorics for the Feynman graphs of 
the anisotropic two-temperature model B is identical to 
the equilibrium model B, and therefore the above renor- 
malization constants assume the same form as their famil- 
iar equilibrium counterparts, apart from the shifted criti- 
cal dimension and a modified geometric factor ^4(d||, d±). 
The one-loop RG flow functions thus are 



Cr 



= Cs = 

n + 2 „ 
■ u 



(> 



B~ = uC~ = u [ — e 



n + ; 



with the stable fixed point [compare Eq. ( 3.72j )] 

5 ^^T8 e + ^' 

leading to the 0(e) exponents, with e = 4 — d — du 

2 



= 2- 



■8 



e + 0(0 



(4.41) 
(4.42) 

(4.43) 



(4.44) 



(4.45) 



For the two- loop fixed point and exponent values to 0(e 2 ) 
for the case n = 1, we refer to Refs. @||. 



4.3 Renormalization of the two-temperature model J 

We now return to the two-temperature model J with non - 
vanis hing mode-coupling term, as defined in Eqs. (1.1C) 
and (4.11). In equilibrium dynamics, the reversible spin 
precession force of model J constitutes a relevant perturba- 
tion to the purely diffusive model B with conserved three- 
component order parameter, and the ensuing non-trivial 
fixed point (with upper critical dimension d^ y — 6) for the 
renormalized effective mode-coupling constant / oc g 2 /A 2 
changes the dynamic critical exponent from z = 4 — r\ 
to z = Jd + 2 — r7)/2, describing faster relaxation pro- 



cesses 



uld expect the reversible mode 
coupling to be relevant in the two-temperature variant of 
model J as well, the issue therefore is, does the RG flow 
again lead to a non-trivial stable fixed point, and what are 
the values of the ensuing critical exponents? 
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It is essential to note, however, that the two-tempera- 
ture model J with go ^ cannot be recast as the dy- 
namics of an eq uivale nt equilibrium model, with an effec- 
tive free e nergy ( 4.13j ), because for the reversible forces in 
Eq. (pOD, 



d d 2 



5S${x) 



0n 




and the necessary Deker-Haake integrability condition 
that would ensure the stability of the equilibrium pro : 
bility distribution oc exp(— H c s[{Sq }]) is violated, except 
for Co = or du = 0. We remark that this is actually a con- 
sequence of the inseparability of "statics" and dynamics 
in the intrinsically dynamic two-temperature model with 
reversible mode couplings. For, if we could first and sep- 
arately consider the effective static free energy, and only 
subs equen tly introduce the dynamics with the analog of 
Eq. (2.12), the elastic pseudo-dipolar propagator would 
appear in the mode-coupling vertex, and the above inte- 
grability condition 



d d 2 



-9o ^ 

0,1 



5H cfi 



-Heff 







(4.47) 

would be satisfied. Following the standard equilibrium pro- 
cedures [ pi| , th e ensu ing cr itical exponents would be given 
by Eqs. ( f4.45| ) and ( 4.36 ), with the exact dynamic ex- 
ponent z = (d + 2 — ry)/2. Yet, such a procedure is not 
possible here, and the derivation of the effective equa- 
tion of motion has t o proceed in the dynamic functional, 
leading to Eq. ( 4.10 ). As opposed to the two- temperature 
model B, the two-temperature model J thus represents 
a genuinely non- equilibrium dynamical model, for, as we 
shall see below, the renormalization-group flow does not 
take the renormalized mode coupling g to zero under scale 
transformations. 

Thus we can invoke no fluctuation-dissipation theo- 
rem in order to relate vertex and response function renor- 
malizations, and we have t o com p ute a lmost all t he Z 
factors, as defined in Eqs. ( |4.15| )-( [ll"7| ) and ( [4.25[ ), in- 
dependently. Neither is there a Ward identity relating 
the renormalization of the mode-coupling vertex to simple 
field renormalizations |p9[ . Fortunately, though, because 
of the momentum dependence of the mode-coupling ver- 

tex ' / os»s^(" qi ' 0;£ 2 L " Pi ' 0; 3 r+P±,Q) = <?o(q_L- 
p ± ) € a ^ + 0(g%), at least Eq. (|4.22|) is valid for the two- 



temperature model J as well. Consequently the identi- 
ties (4.23) still hold, leading immediately to Eq. (4.36) 



or Zu = 2, as to be expected. Furthermore, simple power 
counting yields that the dynamical upper critical dimen- 
sion, where the mode-coupling constant becomes dimen- 
sionless, is 



d^ = 6 - d 



(4.48) 



i.e., the spatial anisotropy reduces d^ y from its equilibrium 
value in exactly the same way by du as the static upper 



critical dimension (4.14). Thus, we define a dimensionless 
renormalized mode-coupling constant according to 

,1/2- 



g = Z^g Q B(d h d ± y^fi- 



-'m--' 2 , (4.49) 
where e = c?J? y — d = 6 — d — du = 6 — 2d» — d±, and 



B(d h d±) = 



r(4-d/2-d\\/2)r(d/2) 



c ^ii/ 2 2^/2 r(d±/2) 



(4.50) 



The appearance of two different upper critical dimensions 
implies that we shall have to compute fixed points and 
exponents in a double expansion with respect to e and e. 
In order to evaluate the renormalization constants, we 



start with the dynamic susceptibility (4.5), and first com- 



pute the fluctuation-induced T c shift from the condition 
Xo(q|| = 0,q^ — * 0,cj = 0) _1 = for r = r 0c - Introduc- 
ing the effective mode-coupling constant 



/o 



38 



2d 



(4.51) 



we find in dimensional regularization 

e - d ~ d ii 5u A(d\\,d ± )\r 0c \ 



roc 



3 {d + d\\-2)(A-d- d\\) 

d\\d± Jo B(d\\,d±_) 

' (d- 2) (4-d-d||)(6-d-d||) ' 



(4.52) 



to be compared with Eq. (4.38) for the two-temperature 
model B. Subsequently rendering both d q 2 xfaw = OjQXj 

^ = 0)- 1 jq i=o and x(q|| = 0,q± -> 0,uj = 0) UV-finite 
yields, after a somewhat tedious calculation, the Z factors 



Z = 1 



rf||(6 - d\\) fg B(dpdj_) ii~ 
2(4 — rf N ) e 
zz =1 5 up A(d\\ , d±) fi~ e 

6 t 
d||(3-d||) JoBid^dx)^ 



(4.53) 



(4.54) 



Here, we have employed minimal subtraction, where only 
the residues of the singular e poles were retained; i.e., in 
the expressions cx /o, d was replaced with 6 — du , and d± 
with 6 — 2d|| to this order. The diffusion constant renor- 
malization is then most conveniently found by considering 



the composite operator <9 g ^ [Ag 2 
0)]| qj _=O) with the result 



■ r s[SS] fall =0,q±,w = 



Z- 1 Z x = 1 



2(2 -d||) f B(d h dx)tt- 



A-d 



(4.55) 



Equivalently, this combination of Z factors can be estab- 
lishes by comparing the UV singularities in d q ^ J~ g (q|| = 
0, qj_,w = 0)| qi= o with the previously established ones 
in x(l|| = 0, q^ — > 0,c<j = 0). From renormalizing the 
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noise vertex function <9 ? 2 r~~(q|| = 0, q±,lu 
we obtain 



Za 1 Z\ 



f B(d\\,d±) [i 



lq_L=0, 



(4.56) 



7df - 76d|| + 48 
6(4 -dii) 



(4.67) 



At last, by means of rather lengthy calculations for the 
derivatives 9( qi . Pl) r-^^-qx, 0; 3^-px,0; 3 j k +px, 
0)| qx= p x=O) 9 ?i r~ sss (-qx,0;^,0;^,Q;Y,0;)| qx=o 
we arrive at the coupling constant renormalizations 



d \\ foB(d\\,d±) fi~ 



Zq Z\ z~ 

tj s\ 



+ 



3 e 
11 u A(d\\ , dj_) /i~ £ 
~ ~6 e 
3-d|| /b -B(d||, dx) A*" 



(4.57) 



(4.58) 



For /o = , Eq s. ( 4.54) an d ( 4.56 ) reduce to the one- loop 
Z factors ( 4.39[ ) and ( 4.40D for the two-temperature model 
B with n = 3, while setting dy — recovers the familiar 
renormalization constants for the equilibrium model J for 
isotropic ferromagnets M.ElU2n] . 



4.4 Discussion of the RG flow equations 

In terms of the renormalized couplings u and 

f = 



2d, A 2 



(4.59) 



the one-loop zeta functions for the two-temperature model 
J become 



Cs 

c = 

Cx- 

(t- 

c~ 



d||(8-d||) 



>s 2(4 -d,|) 

d||(6-d||) ~ 
2(4- d,|) 7 ' 

e d||(16-d||) 
~2 ~ 12(4- d,|) 



'II 



10d,| + 8 



2(4- d,|) 



„ 5 „ 
6 

11 „ 

-e+ — w 
6 



II 



2(4-d||) / ' 
2d||(10- d,|) 
" 3(4 -dy) 



f 



(4.60) 
(4.61) 
(4.62) 
(4.63) 
(4.64) 
(4.65) 



Notice that the dynamic coupling constant enters the RG 
flow for the "static" non-linearity u; this is yet another 
indication that this model is of genuinely dynamical char- 
acter. In the equilibrium limit dii = 0, the statics and 
dynamics decouple. 

The RG beta functions for u and the effective mode 



coupling (4.59) of the two-temperature model J read 



f3 z = u( z = u[-e+—u 



11 „ 2dii(10-d 



3(4 -d,|) 



/ (4.66) 



For d|| = 0, we thus recover the stable equilibrium fixed 
points 



r = § + (e 2 ) 



6 / 9, 

n e+(e) < 



and critical exponents, see Eqs. ( 4.34 ) and ( |4.35| ) 



f] = + O(e 2 ) 



L = 2- 
d+2 



f 0(e 2 ) . 



(4.68) 

(4.69) 
(4.70) 



As can be seen here, the correction to z is merely given 
by the 0(e 2 ) contribution to the static exponent r\ of the 
three-component Heisenberg model. For, in equilibrium 
there is an additional identity Z g = Zg p9f , or £* = 
— (e + rj)/2. The condition for the existence of a non-trivial 
finite fixed point thus becomes C\ — Q = (d — Q ~ ? 7)/2, 
or 

• - ^ ■ (4.TD 



In the full non-equilibrium theory (d\\ > 0), Eq. ( |167| ) 
still implies that the two-temperature model B fixed point 
with /* = is unstable for e = 6 — d— dii > 0. To one- loop 
order, the finite positive fixed point 



/* 



6(4- d 



7df, - 76d,| + 48 



e + 0{e 2 ,e 2 ) 



(4.72) 



exists only in the interval < d|| < |(19— a/277) « 0.6733; 
already for dy = 1, the RG flow takes the mode cou- 
pling to infinity! (According to the one-loop beta func- 
tion (3y, there is another regime where /* > 0, namely 



4 < d|| < f(19 -I 
mensions d? = 4 



v/277) 10.18; yet in suc h hig h di- 
dn < 0.) According to Eq. Q4.66| ), the 





divergence of (3j(£) under scale transformations as £ 

furthermore drives the "static" non-linearity to +oo as 
well. Apparently, the two-temperature model J asymptot- 
ically enters a genuine strong- coupling regime, which does 
not allow for a perturbational calculation of the critical 
exponents (for d\\ > 0.6733) |o). 

Formally though, we may expand about the equilib- 
rium model J, and thus obtain critical exponents in the 
limit d|| <C 1. To first order in due, we find 



2 1+ 3 d > 



6 5 7 

— e + — d\\£ , 



11 



11 



leading to the critical exponents 
3 



n 



■d\\e , 
5 

TT 



Vs 



25 A 

™ due 
66 11 



(4.73) 

(4.74) 
(4.75) 



Uwe C. Tauber et al.: Non-equilibrium critical behavior of 0(n)-symmetric systems 



19 



A £ 1 j 

z = 4 d»e 

2 6 11 

A = l---—d«e. 
4 12 11 



(4.76) 
(4.77) 



Notice, however, that this procedure amounts to an 
expansion with respect to three dimensional parameters, 
namely e = 4 — d — d\\, e = 6 — d — dn, and due. More- 
over, the divergence of the non-expanded fixed point /* at 
dy « 0. 6733 indica tes th at an extrapolation of the formal 
results ( 4.75 ) and ( 4.76 ) to any physical dimension d\\ > 1 
is unlikely to work. On the other hand, we cannot exclude 
that this divergence is merely a one- loop artifact, and is 
cured if one calculates the RG beta functions to higher 
loop orders. Yet another possibility might well be that the 
divergence of /* and u* indicates the absence of a non- 
equilibrium stationary state of the two-temperature model 
J in the vicinity of its critical point. A somewhat less dras- 
tic implication may be that merely perturbation theory 
breaks down, and non-perturbative approaches could pos- 
sibly characterize the scaling behavior at the transition of 
the two-temperature model J successfully. 

At any rate, though, we may draw the following con- 
clusions: (i) As opposed to our non-equilibrium version of 
the SSS model, model J with dynamical anisotropy re- 
mains a genuinely dynamical system, and is very unlikely 
to be described by any simple effective equilibrium the- 
ory at the critical point, certainly not by the equilibrium 
model J. (ii) Obviously the reversible mode coupling term 
in the Langevin equation is highly relevant, driving the 
system away from the well-defined two-temperature model 
B fixed point towards a strong-coupling regime, which at 
least to one-loop order cannot be addressed by means of 
perturbation theory. 



5 Summary and final remarks 

In summary, we have extended previous studies on the 
universality classes of non-equilibrium phase transitions 
by investigating the effect of violating the detailed-balance 
condition in the diffusive dynamics of a conserved field 
which is coupled to the order parameter through reversible 
mode-coupling dynamics. 

In a previous work 0], it was established that the 
universality class of the second-order phase transition in 
a system with non-conserved order parameter dynamics 
is not affected by isotropic breaking of detailed balance. 
Extending this result, we found here that (1) reversible 
mode coupling apparently remains ineffective in generat- 
ing a new universality class even if dynamical anisotropy 
is present in the diffusive dynamics of the conserved field, 
provided the order parameter itself is non-conserved and, 
(2) dynamical anisotropy does become relevant if the or- 
der parameter itself is conserved, in which case reversible 
mode couplings may drive the system towards entirely dif- 
ferent critical behavior, which apparently cannot be de- 
scribed by known equilibrium universality classes. 

These results give further support to previous observa- 
tions that the dynamics of a non-conserved order parame- 



ter (with model A being the simplest realization) is robust 
against non-equilibrium perturbations, while the dynam- 
ics of a conserved order parameter field (model B, if there 
are no reversible mode couplings present) is extremely 
sensitive to detailed-balance breaking through dynamical 
anisotropy in the system. One of the surprising features 
of our results is that dynamical anisotropy, even in com- 
bination with reversible mode coupling terms, cannot de- 
stroy the stability of the equilibrium critical fixed point 
of the non-equilibrium SSS model with non-conserved or- 
der parameter. Only in very extreme cases of effectively 
zero and infinite conserved noise temperatures can one 
find new universality classes, which, however, should well 
influence the crossover behavior and the corrections to 
scaling at the critical point. In our non-equilibrium model 
J, the reversible mode coupling terms had a more dras- 
tic effect, however. While the two-temperature model B 
steady state dynamics can be written in the form of an 
effective equilibrium model with elastic or pseudo-dipolar 
long-range interactions, such a simple representation is not 
possible for the two-temperature model J with dynamical 
anisotropy. We were not even able to identify a stable and 
finite renormalization-group fixed point (to one-loop or- 
der) for this model, but were led to RG runaway flows 
towards a genuinely strong-coupling regime instead (ex- 
cept formally for dn <C 1). 

The fact that the rather complicated combination of 
relaxation, diffusion, reversible mode coupling and dynam- 
ical anisotropy does not change the nature of the second- 
order phase transitions in systems with non-conserved or- 
der parameters, gives us hope that such non-equilibrium 
phase transitions can be understood in terms of a rela- 
tively small number of universality classes. However, for 
systems with conserved order parameters, the situation is 
obviously quite different since both dynamical anisotropy 
itself, and the combination of dynamical anisotropy and 
reversible mode couplings can lead to new universality 
classes. It is clear that the mapping out of the relevant 
non-equilibrium perturbations for the case of conserved 
order parameter remains an open task. It should also be 
noted that most of the studies of this and related problems 
have as yet been restricted to systems with local currents 
(heat baths of different temperatures are attached to the 
system at every point). Problems that are associated with 
global currents, such as, e.g., the driven lattice gas [Q, 
are much richer and more difficult to analyze. Since the 
effective interactions observed here carry dipole-like an- 
gular dependences as well, it may, however, be possible 
that global currents will not cause a significant increase in 
the number of non-equilibrium universality classes. A sys- 
tematic investigation of the effects of global currents, in- 
volving the generalization of the models with local-current 
into global-current models, should clearly be the subject 
of further studies. 
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A Ward identities for the non-equilibrium 
SSS model 

In this appendix we d erive the basic Ward identities that 
were used in Sec. 



3.1 



In the equilibrium SSS model, one 
can use the ordinary Ward identities between the non- 
linear susceptibility and the linear susceptibility jlTj] to 
reduce the number of Z factors required to renormalize 
the theory. However, this procedure is only valid if the re- 
sponse functions renormalize multiplicativcly, a property 
which follows from the fact that in equilibrium the zero- 
frequency limit of a multi-linear response function is equal 
to the corresponding static correlation function [ fl6| . Since 
the non-equilibrium SSS model does not obey detailed bal- 
ance, these Ward identities are not useful in this case. In 
this appendix we will thus work directly with the vertex 
functions, and derive the basic Ward identities which fol- 
low from the 0(n) symmetry of the dynamic functional, 

and the fact that the Mq fields are the generators of this 
symmetry group. 

We consider the following canonical transformation for 
the fields S$ and M Q 



a/3 



6S£ = eA^{Mg' 



5Mf = tA^{M^ ,M% P } 



So) 



e(A av Mf - A^M™) , 



(A.l) 



(A.2) 



where e is a small parameter and A^ v is an arbitrary anti- 
symmetric tensor which is constant in space and time. (In 
this Appendix, we use Einstein's convention of summation 
over repeated indices.) This transformation preserves the 



Poisson brackets between Sq and M 



7<5 



and between M, 



afS 



and Mq s . If this transformation is supplemented with the 

transformation laws for the auxiliary fields Sq and Mq^ , 
as given by 



6S% = eA av S% 

r af3 



SM* P = t{A av Ml 



V0 



#"M m ) 



(A.3) 
(A.4) 



then one can show that the Janssen-de Dominicis func- 
tional J[{Sg},{S$}, {M Q/3 }, UC^I is invariant with re- 
spect to the joint transformation ( |A.l| ) to (A.4). 

If the tensor A^ v (t) is now allowed to depend on time, 
then the dynamic fu nctio nal is no longer invariant under 
the transformations (A^l) to ( A.4 ), but picks up the extra 
terms 



SJ 



d d x 



dtA ap {t) [SS] ap (x,t) 



a/3, 



d d x / dtA af3 {t) [MM] a/3 (x,i) , (A.5) 



where [55]q /3 (x, t) and [MM]q^(x, t) are composite oper- 
ators which are defined as 



1 



[SS] a "(x,t) = 5 [S?(x,t)Sff(x,t)- 
-^(x,f)^(x,i)] 



(A.6) 



and 



[MM]^(x,i) = \ [C(x,t)M^(x,t) - 

-M£"(x,t)Mr(x,t)], (A.7) 

and the "•" stands for differentiation with respect to time. 
In order to proceed, we consider the generating functional 

Z[h a , h a , H a(3 , H a( \ J af3 , J a/3 , C a ] = 

V[{Sg}]V[{^}] { V[{Mf}]V[{iMf}] x 



where we have introduced a source term for the fields Sq , 
S$, Mq 13 , and Mq 13 , [SS]^ and [MM]^ in Z, and also 
a source term for the composite field [MS]g = M a "SJ, 
since this composite operat or al so enters in the definition 
of the linear susceptibility ( |3 .25 ). One can now obtain the 
transformation l aws f or a ny co mposite operator from the 
transformations ( A.l ) to (|A.4| ) 



If one now applies the transformation (A.l) to (A.4) 
with a time-dependent parameter A^ v (t) to the dynamic 
functional Z ', then it is easy to see that the two term s 
generated by the transformation, as given in Eq. (A.5), 
can be absorbed in the transformation law for the sources 
of [SS] f and [MM]f. 

One ends up with the following identity 

Z[h a , h a , H a/3 7 H a/3 7 J a/3 , J af3 7 C a ] = 
Z[h a - eA^h , h a - eA^h 13 , 

H af} ~ e{A av H vfi - A^ v H va ), 

H af3 - e(A au H^ - A^H va ), 

ja (s _ t ^v r s3 _ A P»j»«} _ 2eA af3 , 
C a -eA aP £P] . (A.9) 

Expanding this identity to first order in e, and using 
the fact that A a P(t) is antisymmetric but otherwise ar- 
bitrary, one obtains the following relation for the vertex 
functions 
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= 



(A.10) 



where «g(x,t) = (Sff(x,f)>, Sg>(x,i) = (S Q (x,i)), and 



(x,t) = (M ap (x,i)), m^(x,t) = (M ^(x,i)>. 
A similar identity can be obtained for the correlation 
functions. One can then derive the usual Ward identities 
for multi-linear response functions from this identity. This 
establishes the equivalence of the two procedures in the 
equilibrium case. 



Taking the variational derivative of ( A. 10) with respect 
to m 7l5 (x, t) and m"^(x,(), and setting the source terms 
to zero, one obtains, after taking the Fourier transform, 
the identity 



{' 



0MM[SS] 



(q,w;0,o>) 



0MM[MM] 



(q,w;0,w)J 



(A.ll) 



where we have used the tensor properties of r^ff, ■ This 

. . 0M M . 

identity relates the vertex functions with one insertion of 
the composite operators {SS] a P and [MM] a P to vertex 
functions wi th no insertions. Taking the variational deriva- 
tives of Eq. (A. 10) with respect to the other fields, one can 
obtain similar Ward identities for other vertex functions. 

These Ward identities show that no multiplicative re- 
normalization is needed for the composite operators [5S'] Q ' 3 



[MM] 



1. However, these 



and [MM] af:l , i.e. Z [§s] = Z { 

identities do not exclude the need of additive renormaliza- 
tion for these operators provided that this renormalization 
is oc q 2 

We now consider Eq. (3.25), which follows from the 
identity 



5J 



H af3 (x,t) = , 



5Mf '(M), 



(A.12) 



which in turn can be proven using the fact that the path 
integral of a functional derivative vanishes. If one now 



takes the functional derivative of this equation with re- 
spect to m 7l5 (x', t') a nd p erforms a Fourier transforma- 
tion, one obtains Eq. ( 3.2 5| ) . 

The renormalized version of ( [3.25| ) can be obtained if 
one uses the identity Zjggi = 1 and the definition of the 
renormalized diffusion constant d q 2r^j M (quj) \np= D in 



Eq. ( |3.25D . We thus find 



r MM (q> w ) = iuj + D i 2 + 

+ 2g^A- 1/2 (Z M Z g )- 1 / 2 r Al[ s S] (q,u;) , (A.13) 

where the renormalized vertex function [SS] (l> w ) ^ s 
fined by 



r M[ss](^ w ) - - M 

X ( r 0Af[SS](q^) - 1 2d q 2r 0M[SS]( ( l' UJ ) \np) , 



-1/2 



(A. 14) 



and where the second term on the right-hand-side follows 
from the definition of the renormalized diffusion constant 
and corresponds to an additive renormalization needed to 
render [SS] al3 finite. 

Hence, since all the quantities in ( A.13] ) are finite, we 
conclude that 

Z g Z M = 1 (A. 15) 

must hold to all order s in pe rturb ation theory. If one takes 
the ratio of Eqs. (3.14) and (3.1C), one can check explicitly 
that this identity holds to one loop-order. 

A similar set of Ward identities can be derived for the 
equilibrium model J |^l|,^9) , but it is unclear if these iden- 
tities still hold in the non-equilibrium effective mo del J 
with long-range elastic forces considered in Sec. 4.3. 
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